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Abstract. We investigate classical thermal activation (TA) and macroscopic quantum 
tunneling (MQT) for a Josephson junction coupled with an LC circuit theoretically. The TA 
and MQT escape rate are calculated analytically by taking into account the two-dimensional 
nature of the classical and quantum phase dynamics. We find that the MQT escape rate is 
largely suppressed by the coupling to the LC circuit. On the other hand, this coupling gives 
rise to slight reduction of the TA escape rate. These results are relevant for the interpretation 
of a recent experiment on the MQT and TA phenomena in grain boundary YBCO Josephson 
junctions. 



1. Introduction 

Recent experimental observations of macroscopic quantum tunneling (MQT) [H [H O [H [5] and 
energy level quantization [61 [7] in high-T c superconductor Josephson junctions (JJs) open up 
the possibility for realizing high-T c quantum bits. 

In YBCO grain boundary junctions which was used in MQT experiments [HE], it was found 
that the stray capacitance Cs of the electrodes, due to the large dielectric constant of the STO 
substrate at low temperature (e r > 10000), and the inductance L$, due to the large London 
penetration depth in c-axis and/or Josephson coupling between Cu02 planes in one of the 
electrodes have large influence on the macroscopic dynamics and can be taken into account 
by an extended circuit model (see Fig. 1(a)). In Fig. 1, is the phase difference across the 
Josephson junction and (ft s = (2tt/$>o)IsLs + (ft is the- phase difference across the capacitor Cs, 
where Is is the current through the inductor. It was found that, in the microwave-assisted 
MQT experiment [6], the bias current J ex t dependence of the Josephson plasma frequency u> p is 
quantitatively explained by this model 019]. However, the validity of the extended circuit model 
for TA and MQT escape processes have not yet been explored. Below we calculate the MQT 
and TA escape rate base on the extended circuit model and try to compare with experimental 




results. 



2. Model 

The Hamiltonian of the circuit (Fig. 1(a)) can be written as 



n 



q2j - ^+mo s 



Qi 



2C S 



COS i 



■ + 



{± - <t>sf 
2rj 



i4>s 



(2) 



where Qj = Cj(&o/2ir)(d(j)/dt), Qs = Cs{^o/2ir){d(j)s /dt), Ej is the Josephson coupling 
energy, 77 = 2TrI c L s /®o = L s /Lj (Lj = $o/2vr/c), and 7 = I e xt/^c(4xt is the external 
current). This Hamiltonian describes the quantum dynamics of a Active particle moving in 
the 2-dimensional tilted washboard potential U((j),4>s) (see Fig. 1(b)). By introducing the new 



coordinate (x, y) 



(f) m , 4>s — 4*™) and assuming 7 ~ 1, we can rewrite the Hamiltonian as 
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where 



is the local minimum point, M = Cj(&o/2n) 2 , and m = Cs(^o/2ir)^ 



By using the functional integral method [TOl QTJ Q2J [131 EH 113 US] , the partition function of 
the system can be written as Z = J Vx(t) J 2?y(r)exp — J Q ^ C[x, y]/h , where C[x,y] is the 
Euclidean Lagrangian. The Lagrangian is a quadratic function of y. Therefore the functional 
integral over variable y can be performed explicitly. Then the partition function is reduced to a 
single functional integral over x, i.e., Z = J Vx(t) exp(— S c g[x]/h), where the effective action is 
given by 
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In this equation, Urn = Ej — x 3 /6 + — r f)/2x 
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where lolc = l/v^sCs is the resonance frequency of the LC circuit. Thus the dynamics of 
the phase difference in the 2D potential U((j), 4>s) can be mapped into simple ID model. Note 
that, due to the coupling between the Josephson junction and the LC circuit, S e g [x] contains a 
"dissipation" term (the second term) which is nonlocal in imaginary time r. 



3. Thermal Activation Process 

The thermal activation rate well above the crossover temperature is given by |lL)j 
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where the effective trial frequency wr is a solution of uj 2 ^ + wr7(u;r) 
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(8) 
is the 



Matsubara frequency, to p = y/2irlc /QqCj(1 — 7 2 ) 1//4 is the Josephson plasma frequency, and 
7(u> n ) = {C s / C j)\uj n \oj 2 LC / {oj^ + oj 2 lc ) is the damping function. Note that the potential barrier 
height Vo = (2/27)Mu> 2 x 2 (x\ = 3^1 — 7 2 ) is not modified even in the presence of the LC 
circuit. Therefore the coupling to the LC circuit only modifies the prefactor of Tta- The 
effective trial frequency wr and the quantum correction c qm in the prefactor can be calculated 
analytically, i.e., 
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for the non non-adiabatic cases (wlc *^ w p) which is applicable to the MQT experiment in 
YBCO grain-boundary JJs [H [6]. In the case of Ls/Lj 3> 1, the prefactor of Tta coincides 
with the result without dissipative effects [10J. Therefore we can conclude that the influence of 
the coupling to the LC circuit on the thermal activation process is quite weak. So the system 
behaves as simple ID systems well above the crossover temperature. This result is qualitatively 
consistent with the experimental observation [Jj. 
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4. Macroscopic Quantum Tunneling Process 

The MQT escape rate at zero temperature is defined by Tmqt : 

In order to determine Tmqt, we use the bounce techniques. j.iich wc gci ± M q T 

(ajp/27r)\/l20-KB exp (-B) , where B = (Sq[xb] + Sdiss[ x B])/ft is the bounce exponent, that 
is the value of the action 5 e ff evaluated along the bounce trajectory x_b(t) = xisech(u;pr/2). In 
the non-adiabatic case (ujlc "C u p ), we get 
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where the "dissipation" correction term is given by 5M/M = 5Lj/2Ls- In contrast to T^\, 
the coupling to the LC circuit gives rise to an increase of the bounce exponent. Therefore the 
two-dimensional nature has large influence on the MQT escape process. 

In order to check the validity of our model, we try to compare our result with the experimental 
data of the switching current distribution at the low temperature regime pQ. The switching 
current distribution P(rj) is related to the MQT rate T(7,r/) as 
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where v = \dr//dt\ is the sweep rate of the external bias current. From the experiment pQ, it is 
found that the measured value of the full width at half maximum (HMFW) a of P{rf) is 11.9nA 
at the temperature-independent MQT regime. On the other hand, in order to numerically 
calculate a, we need the information of I c , Cj, L$, and C5. The values of I c = 1.4/xA, 
L s = 1.7nH and Cs = 1.6pF were directly determined from the TA and the microwave-assisted 
MQT experiments [TJ [U]. Therefore the only fitting parameter is Cj. From the numerical 
estimation of a, we found that Cj = 0.2 pF gives good agreement with the experimental value 
of a. The obtained value of Cj is consistent with the estimated value Cj ~ O.IC5 ~ 0.16pF 
based on the geometry of the junction [9]. Therefore, we can conclude that the extended circuit 
model can quantitatively explain the MQT experiment in the YBCO grain boundary junctions. 



5. Summary 

We have theoretically studied the TA and MQT escape for the Josephson junction coupled with 
a LC circuit by taking into account the two-dimensional nature of the phase dynamics. Then 
we found that the coupling to the LC circuit gives negligible reduction for TA escape rate. On 
the other hand, we also found that the MQT escape rate is largely reduced due to the coupling 
with the LC circuit. These results are consistent with experimental data of YBCO Josephson 
junctions. 
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